Binding interaction of laccases from Bacillus Subtilis after industrial dyes exposure: Molecular docking and molecular dynamics simulation studies

Dyes are becoming more widely used around the world wide, but there is no effective bioremediation approach for removing them completely from the environment. Several dyes are mentioned to be degraded through bacteria; however, it's still unknown how the particular enzymes act throughout the dye degradation. The behavior and function of these enzymes in the biodegradation of azo dyes (Textile dyes) had been investigated experimentally by the numbers of the researchers, however, the molecular mechanisms remain unclear. Therefore, the interaction mechanisms of textile dye (methyl orange) with laccase from B. subtilis were explored through molecular docking and molecular dynamics simulations, the three selected dyes (methyl orange, malachite green, and acid blue 62) that interact positively with laccase on the basis of their maximum binding energy, molecular docking results indicate that one of the three dyes is more stable as a target for degradation through Bacillus subtilis laccase. Therefore, subsequent research focused solely on one substrate: methyl orange. Molecular Dynamics simulation study was applied after the molecular docking to determine the interaction between laccases and methyl orange dyes. The trajectory was proved with root mean square deviation and root mean square fluctuation analysis. According to the molecular dynamics simulation results, laccase-methyl orange complexes remain stable during the catalytic reaction. So, this study demonstrates how laccase is involved in methyl orange bioremediation.

Dynamics and Energetics of a Bromodomain in complex with bromosporine from Leishmania donovani

Leishmaniasis continues to be a neglected tropical disease, affecting people and animals and causing significant economic losses. Therefore, there is interest in the study and evaluation of new drug targets. In fact, it has been shown that by interfering with lysine-reading proteins such as bromodomain (BMD) there is a decrease in parasite survival. In this study, we researched the dynamics and energetics of the Leishmania donovani BMD in complex with bromosporin, which is considered to be a pan-inhibitor of BMDs, with the aim of understanding the molecular recognition mechanism. Molecular dynamics (MD) and non-equilibrium free energy calculation guided by steered molecular dynamics (SMD) simulations showed that the BMD has three flexible amino acid regions and bromosporin exhibiting various recognition states during the interaction. These results corroborate the promiscuity of bromosporin for energetically favourable sites, with the possibility of expanding its inhibition to other bromodomains. Furthermore, these results suggest that Van der Waals interactions have more relevance for complex recognition and residues ASN-87 and TRP-93 are key in forming hydrophobic and H-bond interactions, respectively. This research provides new insights for understanding the recognition mechanism, dynamics and energetics of the complex for the development of new therapeutic strategies.

La leishmaniasis sigue siendo una enfermedad tropical desatendida, que afecta a personas y animales y causa importantes pérdidas económicas. De ahí el interés por estudiar y evaluar nuevas dianas farmacológicas. De hecho, se ha demostrado que al interferir con proteínas lectoras de lisina como el bromodominio ("bromodomain", BMD) se produce una disminución de la supervivencia del parásito. En este artículo estudiamos la dinámica y la energética del BMD de Leishmania donovani en complejo con bromosporina, que se considera un pan-inhibidor de BMD, con el objetivo de comprender el mecanismo de reconocimiento molecular. Las simulaciones de dinámica molecular (DM) y el cálculo de energía libre de no-equilibrio guiado por dinámica molecular de estiramiento (DMS) mostraron que BMD tiene tres regiones de aminoácidos flexibles y la bromosporina presenta varios estados de reconocimiento durante la interacción. Estos resultados corroboran la promiscuidad de la bromosporina por sitios energéticamente favorables, siendo posible expandir su inhibición a otros bromodominios. Además, los resultados sugieren que las interacciones de Van der Waals tienen más relevancia para el reconocimiento del complejo y los residuos ASN-87 y TRP-93 son clave en la formación de interacciones hidrofóbicas y de puentes de hidrógeno, respectivamente. Esta investigación proporciona nuevos conocimientos para comprender el mecanismo de reconocimiento molecular, la dinámica y la energética del complejo para el desarrollo de nuevas estrategias terapéuticas.

A leishmaniose continua a ser uma doença tropical negligenciada, que afeta os seres humanos e os animais e causa perdas econômicas significativas. Daí o interesse em estudar e avaliar novos alvos de medicamentos. De fato, a interferência com proteínas leitoras de lisina, como o bromo domínio ("bromodomain", BMD), tem demonstrado diminuir a sobrevivência do parasita. Neste trabalho, estudamos a dinâmica e a energética do BMD de Leish-mania donovani em complexo com a bro-mosporina, considerada um pan-inibidor da BMDs, com o objetivo de compreender o mecanismo de reconhecimento molecular. As simulações de dinâmica molecular (MD) e cálculo de energia livre de não-equilíbrio guiada por dinâmica molecular esticamento (MDS) mostraram que o BMD tem três regiões de aminoácidos flexíveis e que a bromosporina apresenta vários estados de reconhecimento durante a interação. Esses resultados corroboram a promiscuidade da bromosporina para locais energeticamente favoráveis, possibilitando a expansão de sua inibição para outros bromodomínios. Além disso, os resultados sugerem que as interações de Van der Waals são mais relevantes no momento do reconhecimento do complexo e os resíduos ASN-87 e TRP-93 são fundamentais na formação de interações hidrofóbicas e de ligações de hidrogênio, respectivamente. Essa pesquisa fornece novos insights para compreender o mecanismo de reconhecimento, a dinâmica e a energética do complexo para o desenvolvimento de novas estratégias terapêuticas.

A threefold approach including quantum chemical, molecular docking and molecular dynamic studies to explore the natural compounds from Centaurea jacea as the potential inhibitors for COVID-19 / Uma abordagem tríplice incluindo química quântica, docking molecular e estudos dinâmicos moleculares para explorar os compostos naturais de Centaurea jacea como inibidores potenciais para Covid-19

Abstract In the current report, we studied the possible inhibitors of COVID-19 from bioactive constituents of Centaurea jacea using a threefold approach consisting of quantum chemical, molecular docking and molecular dynamic techniques. Centaurea jacea is a perennial herb often used in folk medicines of dermatological complaints and fever. Moreover, anticancer, antioxidant, antibacterial and antiviral properties of its bioactive compounds are also reported. The Mpro (Main proteases) was docked with different compounds of Centaurea jacea through molecular docking. All the studied compounds including apigenin, axillarin, Centaureidin, Cirsiliol, Eupatorin and Isokaempferide, show suitable binding affinities to the binding site of SARS-CoV-2 main protease with their binding energies -6.7 kcal/mol, -7.4 kcal/mol, -7.0 kcal/mol, -5.8 kcal/mol, -6.2 kcal/mol and -6.8 kcal/mol, respectively. Among all studied compounds, axillarin was found to have maximum inhibitor efficiency followed by Centaureidin, Isokaempferide, Apigenin, Eupatorin and Cirsiliol. Our results suggested that axillarin binds with the most crucial catalytic residues CYS145 and HIS41 of the Mpro, moreover axillarin shows 5 hydrogen bond interactions and 5 hydrophobic interactions with various residues of Mpro. Furthermore, the molecular dynamic calculations over 60 ns (6×106 femtosecond) time scale also shown significant insights into the binding effects of axillarin with Mpro of SARS-CoV-2 by imitating protein like aqueous environment. From molecular dynamic calculations, the RMSD and RMSF computations indicate the stability and dynamics of the best docked complex in aqueous environment. The ADME properties and toxicity prediction analysis of axillarin also recommended it as safe drug candidate. Further, in vivo and in vitro investigations are essential to ensure the anti SARS-CoV-2 activity of all bioactive compounds particularly axillarin to encourage preventive use of Centaurea jacea against COVID-19 infections.

Resumo No presente relatório, estudamos os possíveis inibidores de Covid-19 de constituintes bioativos de Centaurea jacea usando uma abordagem tripla que consiste em técnicas de química quântica, docking molecular e dinâmica molecular. Centaurea jacea é uma erva perene frequentemente usada em remédios populares de doenças dermatológicas e febre. Além disso, as propriedades anticâncer, antioxidante, antibacteriana e antiviral de seus compostos bioativos também são relatadas. A Mpro (proteases principais) foi acoplada a diferentes compostos de Centaurea jacea por meio de docking molecular. Todos os compostos estudados, incluindo apigenina, axilarina, Centaureidina, Cirsiliol, Eupatorina e Isokaempferide, mostram afinidades de ligação adequadas ao sítio de ligação da protease principal SARS-CoV-2 com suas energias de ligação -6,7 kcal / mol, -7,4 kcal / mol, - 7,0 kcal / mol, -5,8 kcal / mol, -6,2 kcal / mol e -6,8 kcal / mol, respectivamente. Dentre todos os compostos estudados, a axilarina apresentou eficiência máxima de inibidor, seguida pela Centaureidina, Isokaempferida, Apigenina, Eupatorina e Cirsiliol. Nossos resultados sugeriram que a axilarina se liga aos resíduos catalíticos mais cruciais CYS145 e HIS41 do Mpro, além disso a axilarina mostra 5 interações de ligações de hidrogênio e 5 interações hidrofóbicas com vários resíduos de Mpro. Além disso, os cálculos de dinâmica molecular em uma escala de tempo de 60 ns (6 × 106 femtossegundos) também mostraram percepções significativas sobre os efeitos de ligação da axilarina com Mpro de SARS-CoV-2 por imitação de proteínas como o ambiente aquoso. A partir de cálculos de dinâmica molecular, os cálculos RMSD e RMSF indicam a estabilidade e dinâmica do melhor complexo ancorado em ambiente aquoso. As propriedades ADME e a análise de previsão de toxicidade da axilarina também a recomendaram como um candidato a medicamento seguro. Além disso, as investigações in vivo e in vitro são essenciais para garantir a atividade anti-SARS-CoV-2 de todos os compostos bioativos, particularmente a axilarina, para encorajar o uso preventivo de Centaurea jacea contra infecções por Covid-19.

Advances in applications of molecular dynamics simulation techniques to the research of self-assembled nano-drug delivery systems / 药学学报

Molecular dynamics simulation technology relies on Newtonian mechanics to simulate the motion of molecular system of the real system by computer simulation. It has been used in the research of self-assembly processes illustration and macroscopic performance prediction of self-assembly nano-drug delivery systems (NDDS) in recent years, which contributes to the facilitation and accurate design of preparations. In this review, the definitions, catalogues, and the modules of molecular dynamics simulation techniques are introduced, and the current status of their applications are summarized in the acquisition and analysis of microscale information, such as particle size, morphology, the formation of microdomains, and molecule distribution of the self-assembly NDDS and the prediction of their macroscale performances, including stability, drug loading capacity, drug release kinetics and transmembrane properties. Moreover, the existing applications of the molecular dynamic simulation technology in the formulation prediction of self-assembled NDDS were also summarized. It is expected that the new strategies will promote the prediction of NDDS formulation and lay a theoretical foundation for an appropriate approach in NDDS studies and a reference for the wider application of molecular dynamics simulation technology in pharmaceutics.

Fosfatase Cdc25B: Interações transientes intramoleculares e complexação com pequenos ligantes / Cdc25B phosphatase: Transient intramolecular interactions and small molecules complexation

Proteínas tirosina-fosfatase (PTPs) possuem papel fundamental na regulação da transdução de sinais e estão envolvidas em diversos processos fundamentais do ciclo celular. As Cdc25 (Cell Division Cycle 25) são fosfatases duais encontradas em todos os organismos eucarióticos e atuam em checkpoints do ciclo celular, permitindo ou inibindo o prosseguimento deste. Este grupo de proteínas pertence à classe de PTPs com atividade baseada em cisteína, apresenta domínio catalítico altamente conservado assim como o motivo catalítico, P-loop. Devido sua função, as Cdc25 são consideradas possíveis alvos terapêuticos para tratamento de câncer e sua interação com pequenas moléculas e inibidores tem sido investigada de forma que análises estruturais e de ligação das Cdc25 com inibidores podem elucidar aspectos importantes do mecanismo de ação destes além de direcionar para o desenho racional de fármacos. Interações cátion-π são interações intra ou intermoleculares não-covalentes que ocorrem entre uma espécie química catiônica, como o grupo guanidino de argininas, e uma das faces de um sistema π rico em elétrons, como dos anéis indólicos de triptofanos. Apesar de pouco discutidas na literatura, quando em comparação às interações não-covalentes mais convencionais, do ponto de vista energético as interações cátion-π são tão importantes na estruturação de proteínas quanto às ligações de hidrogênio ou pontes salinas. De fato estas interações são observadas com frequência em estruturas proteicas resolvidas. O domínio catalítico da Cdc25B possui diversas argininas expostas em sua superfície e um único resíduo de triptofano localizado na região C-terminal flexível, muito próximo do sítio catalítico da proteína. A flexibilidade de proteínas ou de regiões proteicas apresenta importante papel no reconhecimento entre biomoléculas participantes de vias de sinalização e tem sido muito estudada atualmente. Aqui, simulações de dinâmica molecular, experimentos de 1H-15N HSQC RMN, ensaios de cinética de inibição e de ancoragem molecular, evidenciam a existência de contatos cátion-π transientes na superfície de um importante membro da família das Cdc25, a Cdc25B, e de sítios de interação entre inibidores testados e a proteína com destaque a sítios na proximidades do P-loop, região próxima ao C-terminal desordenado, onde se demonstra estabilidade da interação com os pequenos ligantes

Protein tyrosine phosphatase (PTPs) play a fundamental role in the regulation of signal transduction and are involved in several fundamental processes of the cell cycle. Cdc25 (Cell Division Cycle 25) are dual phosphatases found in all eukaryotic organisms and act at checkpoints of the cell cycle, allowing or inhibiting its progression. This group of proteins belongs to the class of PTPs with cysteine-based activity, presenting a highly conserved catalytic domain as well as the catalytic motif, P-loop. Due to their function, Cdc25 are considered possible therapeutic targets for cancer treatment and their interaction with small molecules and inhibitors has been investigated so that structural and binding analyzes of Cdc25 with inhibitors can elucidate important aspects of their mechanism of action besides directing to rational drug design. Cation-π interactions are non-covalent intra- or intermolecular interactions that occur between a cationic chemical species, such as the guanidino group of arginines, and one of the faces of an electron-rich system, such as the indole rings of tryptophans. Although little discussed in the literature, when compared to more conventional non-covalent interactions, from the energetic point of view, cation-π interactions are as important in the structuring of proteins as hydrogen bonds or salt bridges. In fact, these interactions are frequently observed in solved protein structures. The catalytic domain of Cdc25B has several arginines exposed on its surface and a single tryptophan residue located in the flexible C-terminal region, very close to the catalytic site of the protein. The flexibility of proteins or protein regions plays an important role in the recognition between biomolecules participating in signaling pathways and has been extensively studied today. Here, molecular dynamics simulations, 1H-15N HSQC NMR experiments, inhibition kinetics and molecular anchoring assays, evidence the existence of transient cation-π contacts on the surface of an important member of the Cdc25 family, Cdc25B, and of sites of interaction between tested inhibitors and the protein, with emphasis on sites in the vicinity of the P-loop, a region close to the disordered C-terminus, where stability of the interaction with the small ligands is demonstrated

Mechanical properties and immunotherapeutic effects of dissolving microneedles with different drug loadings based on hyaluronic acid

Abstract Improving vaccine immunity and reducing antigen usage are major challenges in the clinical application of vaccines. Microneedles have been proven to be painless, minimally invasive, highly efficient, and have good patient compliance. Compared with traditional transdermal drug delivery, it can effectively deliver a large-molecular-weight drug into the skin, resulting in a corresponding immune response. However, few studies have examined the relationship between microneedle loading dose and immune effects. In this study, the hyaluronic acid (HA) conical and pyramidal dissolving microneedles were prepared by the two-step vacuum drying method, respectively. The model drug ovalbumin (OVA) was added to HA to prepare dissolving microneedles with different loading amounts. The mass ratios of HA to OVA were 5:1, 5:3, and 5:5. The mechanical properties of the dissolving microneedles were characterized using nanoindentation and in vitro puncture studies. The immune effects of the matrix and drug content were studied in Sprague-Dawley (SD) rats. Finally, the diffusion behavior of OVA and the binding mode of HA and OVA in the microneedles were simulated using Materials Studio and Autodocking software. The experimental results showed that the conical microneedles exhibited better mechanical properties. When the mass ratio of HA to OVA was 5:3, the immune effect can be improved by 37.01% compared to subcutaneous injection, and achieved a better immune effect with relatively fewer drugs. This conclusion is consistent with molecular simulations. This study provides theoretical and experimental support for the drug loading and efficacy of microneedles with different drug loadings

A threefold approach including quantum chemical, molecular docking and molecular dynamic studies to explore the natural compounds from Centaurea jacea as the potential inhibitors for COVID-19 / Uma abordagem tríplice incluindo química quântica, docking molecular e estudos dinâmicos moleculares para explorar os compostos naturais de Centaurea jacea como inibidores potenciais para Covid-19

In the current report, we studied the possible inhibitors of COVID-19 from bioactive constituents of Centaurea jacea using a threefold approach consisting of quantum chemical, molecular docking and molecular dynamic techniques. Centaurea jacea is a perennial herb often used in folk medicines of dermatological complaints and fever. Moreover, anticancer, antioxidant, antibacterial and antiviral properties of its bioactive compounds are also reported. The Mpro (Main proteases) was docked with different compounds of Centaurea jacea through molecular docking. All the studied compounds including apigenin, axillarin, Centaureidin, Cirsiliol, Eupatorin and Isokaempferide, show suitable binding affinities to the binding site of SARS-CoV-2 main protease with their binding energies -6.7 kcal/mol, -7.4 kcal/mol, -7.0 kcal/mol, -5.8 kcal/mol, -6.2 kcal/mol and -6.8 kcal/mol, respectively. Among all studied compounds, axillarin was found to have maximum inhibitor efficiency followed by Centaureidin, Isokaempferide, Apigenin, Eupatorin and Cirsiliol. Our results suggested that axillarin binds with the most crucial catalytic residues CYS145 and HIS41 of the Mpro, moreover axillarin shows 5 hydrogen bond interactions and 5 hydrophobic interactions with various residues of Mpro. Furthermore, the molecular dynamic calculations over 60 ns (6×106 femtosecond) time scale also shown significant insights into the binding effects of axillarin with Mpro of SARS-CoV-2 by imitating protein like aqueous environment. From molecular dynamic calculations, the RMSD and RMSF computations indicate the stability and dynamics of the best docked complex in aqueous environment. The ADME properties and toxicity prediction analysis of axillarin also recommended it as safe drug candidate. Further, in vivo and in [...].

No presente relatório, estudamos os possíveis inibidores de Covid-19 de constituintes bioativos de Centaurea jacea usando uma abordagem tripla que consiste em técnicas de química quântica, docking molecular e dinâmica molecular. Centaurea jacea é uma erva perene frequentemente usada em remédios populares de doenças dermatológicas e febre. Além disso, as propriedades anticâncer, antioxidante, antibacteriana e antiviral de seus compostos bioativos também são relatadas. A Mpro (proteases principais) foi acoplada a diferentes compostos de Centaurea jacea por meio de docking molecular. Todos os compostos estudados, incluindo apigenina, axilarina, Centaureidina, Cirsiliol, Eupatorina e Isokaempferide, mostram afinidades de ligação adequadas ao sítio de ligação da protease principal SARS-CoV-2 com suas energias de ligação -6,7 kcal / mol, -7,4 kcal / mol, - 7,0 kcal / mol, -5,8 kcal / mol, -6,2 kcal / mol e -6,8 kcal / mol, respectivamente. Dentre todos os compostos estudados, a axilarina apresentou eficiência máxima de inibidor, seguida pela Centaureidina, Isokaempferida, Apigenina, Eupatorina e Cirsiliol. Nossos resultados sugeriram que a axilarina se liga aos resíduos catalíticos mais cruciais CYS145 e HIS41 do Mpro, além disso a axilarina mostra 5 interações de ligações de hidrogênio e 5 interações hidrofóbicas com vários resíduos de Mpro. Além disso, os cálculos de dinâmica molecular em uma escala de tempo de 60 ns (6 × 106 femtossegundos) também mostraram percepções significativas sobre os efeitos de ligação da axilarina com Mpro de SARS-CoV-2 por imitação de proteínas como o ambiente aquoso. A partir de cálculos de dinâmica molecular, os cálculos RMSD e RMSF indicam a estabilidade e dinâmica do melhor complexo ancorado em ambiente [...].

A threefold approach including quantum chemical, molecular docking and molecular dynamic studies to explore the natural compounds from Centaurea jacea as the potential inhibitors for COVID-19

Abstract In the current report, we studied the possible inhibitors of COVID-19 from bioactive constituents of Centaurea jacea using a threefold approach consisting of quantum chemical, molecular docking and molecular dynamic techniques. Centaurea jacea is a perennial herb often used in folk medicines of dermatological complaints and fever. Moreover, anticancer, antioxidant, antibacterial and antiviral properties of its bioactive compounds are also reported. The Mpro (Main proteases) was docked with different compounds of Centaurea jacea through molecular docking. All the studied compounds including apigenin, axillarin, Centaureidin, Cirsiliol, Eupatorin and Isokaempferide, show suitable binding affinities to the binding site of SARS-CoV-2 main protease with their binding energies -6.7 kcal/mol, -7.4 kcal/mol, -7.0 kcal/mol, -5.8 kcal/mol, -6.2 kcal/mol and -6.8 kcal/mol, respectively. Among all studied compounds, axillarin was found to have maximum inhibitor efficiency followed by Centaureidin, Isokaempferide, Apigenin, Eupatorin and Cirsiliol. Our results suggested that axillarin binds with the most crucial catalytic residues CYS145 and HIS41 of the Mpro, moreover axillarin shows 5 hydrogen bond interactions and 5 hydrophobic interactions with various residues of Mpro. Furthermore, the molecular dynamic calculations over 60 ns (6×106 femtosecond) time scale also shown significant insights into the binding effects of axillarin with Mpro of SARS-CoV-2 by imitating protein like aqueous environment. From molecular dynamic calculations, the RMSD and RMSF computations indicate the stability and dynamics of the best docked complex in aqueous environment. The ADME properties and toxicity prediction analysis of axillarin also recommended it as safe drug candidate. Further, in vivo and in vitro investigations are essential to ensure the anti SARS-CoV-2 activity of all bioactive compounds particularly axillarin to encourage preventive use of Centaurea jacea against COVID-19 infections.

Resumo No presente relatório, estudamos os possíveis inibidores de Covid-19 de constituintes bioativos de Centaurea jacea usando uma abordagem tripla que consiste em técnicas de química quântica, docking molecular e dinâmica molecular. Centaurea jacea é uma erva perene frequentemente usada em remédios populares de doenças dermatológicas e febre. Além disso, as propriedades anticâncer, antioxidante, antibacteriana e antiviral de seus compostos bioativos também são relatadas. A Mpro (proteases principais) foi acoplada a diferentes compostos de Centaurea jacea por meio de docking molecular. Todos os compostos estudados, incluindo apigenina, axilarina, Centaureidina, Cirsiliol, Eupatorina e Isokaempferide, mostram afinidades de ligação adequadas ao sítio de ligação da protease principal SARS-CoV-2 com suas energias de ligação -6,7 kcal / mol, -7,4 kcal / mol, - 7,0 kcal / mol, -5,8 kcal / mol, -6,2 kcal / mol e -6,8 kcal / mol, respectivamente. Dentre todos os compostos estudados, a axilarina apresentou eficiência máxima de inibidor, seguida pela Centaureidina, Isokaempferida, Apigenina, Eupatorina e Cirsiliol. Nossos resultados sugeriram que a axilarina se liga aos resíduos catalíticos mais cruciais CYS145 e HIS41 do Mpro, além disso a axilarina mostra 5 interações de ligações de hidrogênio e 5 interações hidrofóbicas com vários resíduos de Mpro. Além disso, os cálculos de dinâmica molecular em uma escala de tempo de 60 ns (6 × 106 femtossegundos) também mostraram percepções significativas sobre os efeitos de ligação da axilarina com Mpro de SARS-CoV-2 por imitação de proteínas como o ambiente aquoso. A partir de cálculos de dinâmica molecular, os cálculos RMSD e RMSF indicam a estabilidade e dinâmica do melhor complexo ancorado em ambiente aquoso. As propriedades ADME e a análise de previsão de toxicidade da axilarina também a recomendaram como um candidato a medicamento seguro. Além disso, as investigações in vivo e in vitro são essenciais para garantir a atividade anti-SARS-CoV-2 de todos os compostos bioativos, particularmente a axilarina, para encorajar o uso preventivo de Centaurea jacea contra infecções por Covid-19.

A computational predicting of possible inhibitors of the main SARS-CoV-2 protease found in Algerian herbal medicines

COVID-19 is a zoonotic viral disease caused by the SARS-CoV-2 virus. Its abrupt outbreak has caused a tremendous challenge to public health systems due to the rapid spread of the virus. In this sense, a great deal of work has been focused on finding substances from herbal plants to be used against this virus. In order to investigate the molecular interactions between natural metabolites from Algerian herbal plants and the SARS-CoV-2 protease Mpro, computational docking and molecular dynamics were used, also the drug likeness degree and in silico ADMET prediction were carried out in this study. warfarin and catalponol preferentially binds to a pocket of the SARS-Cov-2 Mpro active site that is made up of residues His 41 to Glu 166 and Leu 27 to His 163 with a relatively low binding energy of -7.1 and -6.6 kcal/mol respectively. Dynamic molecular assay further established that only warfarin managed to stay in the active site. The results suggest that warfarin may be an interesting candidate for development as a medical treatment of COVID-19 and more research is proposed, without disregarding its toxicity which deserves to be well studied.

El COVID-19 es una enfermedad zoonótica causada por el virus SARS-CoV-2. Su abrupto brote en años recientes ha supuesto un tremendo desafío para los sistemas de salud pública, como resultado de la rápida propagación del virus. En tal sentido, muchos trabajos se han centrado en encontrar sustancias de origen vegetal, para ser utilizadas contra este virus. Se realizaron estudios de acoplamiento computacional y dinámica molecular para investigar las interacciones moleculares entre los metabolitos secundarios de las plantas herbales argelinas con la Proteasa Mpro del SARS-CoV-2, también se realizaron estudios de semejanza con drogas mediante ADMET computacional. La warfarina y el catalponol se unen preferentemente al sitio activo SARS-Cov-2 Mpro que se compone de residuos His 41 a Glu 166 y Leu 27 a His 163 con una energía de enlace relativamente baja, -7,1 y -6,6 kcal/mol respectivamente. Los ensayos de dinámica molecular establecieron además que sólo la warfarina logró permanecer en el sitio activo. Estos resultados sugieren que la warfarina puede ser un candidato interesante para el desarrollo como tratamiento médico de COVID-19 e instan a realizar más investigaciones, sin dejar de lado estudios de toxicidad respectivos.

A COVID-19 é uma doença zoonótica causada pelo vírus SARS-CoV-2, cujo surto abrupto nos últimos anos representou um tremendo desafio para os sistemas de saúde pública devido à rápida disseminação do vírus. Nesse sentido, muitos trabalhos têm se concentrado em encontrar substâncias de origem vegetal, para serem utilizadas contra esse vírus. Estudos de ancoragem computacional e dinâmica molecular foram conduzidos para investigar as interações moleculares entre metabólitos secundários de ervas argelinas com o SARS-CoV-2 Protease Mpro, estudos de similaridade de drogas também foram conduzidos usando ADMET in silico. A varfarina e o catalponol ligam-se preferencialmente ao sítio ativo SARS-Cov-2 Mpro que é composto pelos resíduos His 41 a Glu 166 e Leu 27 a His 163 com uma energia de ligação relativamente baixa, -7,1 e -6,6 kcal/mol, respectivamente. Ensaios de dinâmica molecular estabeleceram ainda que apenas a varfarina conseguiu permanecer no sítio ativo. Esses resultados sugerem que a varfarina pode ser um candidato interessante para desenvolvimento como tratamento médico para COVID-19 e exigem mais pesquisas, incluindo os respectivos estudos de toxicidade.

Cryptosporidiosis, a public health challenge: A combined 3D shape-based virtual screening, docking study, and molecular dynamics simulation approach to identify inhibitors with novel scaffolds for the treatment of cryptosporidiosis

Cryptosporidiosis is a neglected tropical disease caused by the protozoan parasite Cryptosporidium parvum. Limited therapeutic options, limitation in in vitro parasite culture, and lack of a reliable animal model of parasite for replication of in vivo life cycle and drug testing demand alternative methods for drug development. The in silico methods of drug discovery prove a crucial process in such conditions.Recent research reported a limited number of small molecules for drug development. Purine nucleotide biosynthesis in Cryptosporidium species is dependent on the IMPDH (CpIMPDH) enzyme, so distortion of parasite IMPDH has been pursued as a compelling strategy for curbing Cryptosporidium infection due to its different kinetics from the host enzyme. Our study's primary aim was to discover novel ligand molecules with noticeable activity against Cryptosporidium parvum IMPDH. For this purpose, we selected 18 previously discovered ligands to understand the interaction feature between ligand and receptor, and their shape and electronic features are employed as a template for shape-based virtual screening of the ZINC database (drug-like subset) search approach via Schrodinger-2019 (Maestro 11.9). The obtained hits were subsequently subjected to structure-based screening, quantum polarized ligand docking (QPLD), and molecular dynamics simulations to fetch potential small molecules with the highest binding affinity for CpIMPDH protein. Further ligand binding energy and pharmacokinetic analysis were also taken into consideration as filtering criteria for selecting the most promising drug-like compounds. On this experimentation analysis, three top-ranked (ZINC24855054, ZINC58171263, and ZINC08000072) molecules were found to have appropriate pharmacokinetic properties along with surpassing in silico inhibitory potential towards the CpIMPDH compared to known inhibitors. The molecular docking and molecular dynamics simulation analysis results satisfactorily confirmed the inhibitory action. Therefore, these new scaffolds deduced by the presented computational methodology could recommend lead molecules for designing promising anti-cryptosporidial drugs targeting CpIMPDH protein.

Functional and binding studies of gallic acid showing platelet aggregation inhibitory effect as a thrombin inhibitor / 中草药·英文版

Objective: This study was devoted to identifying natural thrombin inhibitors from traditional Chinese medicine (TCM) and evaluating its biological activity in vitro and binding characteristics. Methods: A combination strategy containing molecular docking, thrombin inhibition assay, surface plasmon resonance (SPR) and molecular dynamics simulation were applied to verify the study result. Results: Gallic acid was confirmed as a direct thrombin inhibitor with IC

Mechanism of Astragaloside Ⅳ on HepG2 Cells Based on Molecular Dynamics Simulation and Experimental Evaluation / 肿瘤防治研究

Objective To reveal the mechanism of action of AS-Ⅳ on HepG2 cells based on molecular dynamics simulation and experimental evaluation. Methods We constructed a "drug-disease" network pharmacological map, analyzed the core genes of astragaloside Ⅳ (AS-Ⅳ) in HCC, screened key signaling pathways, and established a "drug-target" molecular dynamics model. In vitro assay was used to detect migration, proliferation and invasion abilities. Flow cytometry and qRT-PCR were used to detect the effect of AS-Ⅳ on the cell cycle and apoptosis, and the expression of core gene of HepG2. Results The core target of AS-Ⅳ acting on HCC was VEGFA. Compared with the control group, the high concentration of AS-Ⅳ significantly inhibited the migration, invasion and proliferation of HepG2 cells, blocked the metastasis of HepG2 cells from G1 to G2 phase, promoted their apoptosis, down-regulated VEGFA expression and up-regulated TGF-β1 expression. Conclusion AS-Ⅳ may inhibit the proliferation of hepatocellular carcinoma cells through multi-target and multi-pathway.

Computational approaches towards the discovery and optimisation of cruzain inhibitors

The need to develop safer and more efficacious drugs to treat Chagas disease has motivated the search for cruzain inhibitors. Cruzain is the recombinant, truncated version of cruzipain, a cysteine protease from Trypanosoma cruzi with important roles during the parasite life cycle. Several computational techniques have been applied to discover and optimise cruzain inhibitors, providing a molecular basis to guide this process. Here, we review some of the most recent computational studies that provided important information for the design of cruzain inhibitors. Moreover, we highlight the diversity of applications of in silico techniques and their impact.

Computational analysis of fusion protein of anti-HER2 scFv and alpha luffin: A new immunotoxin protein for HER2 positive cancers

Abstract The present study deals with the computational design and analysis of a novel fusion protein based on a single chain variable fragment that binds to the extracellular domain of human epidermal growth factor receptor 2 (HER2) in breast cancer cells. Alpha luffin, a small ribosome inactivating protein (RIP), was attached to the anti-HER2 antibody fragment. I-TASSER modeling provided the full-length structure of the fusion protein. Molecular docking evaluated the molecular interactions of the complementarity-determining regions of designed fusion protein to HER2. Energy minimization and molecular dynamics simulations were conducted to refine the complexes. RMSD plot revealed reasonable stability of the fusion protein during the simulation. The free binding energy profile of complexes affirmed a favorable binding affinity of proteins in complex with HER2 using molecular mechanics Poisson-Boltzmann surface area (G-MMPBSA) algorithm. In general, this approach looks promising in the development of new fusion proteins in terms of immunotoxins with appropriate cytotoxicity.

Improving the thermal stability of Proteus mirabilis lipase based on multiple computational design strategies / 生物工程学报

Proteus mirabilis lipase (PML) features tolerance to organic solvents and great potential for biodiesel synthesis. However, the thermal stability of the enzyme needs to be improved before it can be used industrially. Various computational design strategies are emerging methods for the modification of enzyme thermal stability. In this paper, the complementary algorithm-based ABACUS, PROSS, and FoldX were employed for positive selection of PML mutations, and their pairwise intersections were further subjected to negative selection by PSSM and GREMLIN to narrow the mutation library. Thereby, 18 potential single-point mutants were screened out. According to experimental verification, 7 mutants had melting temperature (Tm) improved, and the ΔTm of K208G and G206D was the highest, which was 3.75 ℃ and 3.21 ℃, respectively. Five mutants with activity higher than the wild type (WT) were selected for combination by greedy accumulation. Finally, the Tm of the five-point combination mutant M10 increased by 10.63 ℃, and the relative activity was 140% that of the WT. K208G and G206D exhibited certain epistasis during the combination, which made a major contribution to the improvement of the thermal stability of M10. Molecular dynamics simulation indicated that new forces were generated at and around the mutation sites, and the rearrangement of forces near G206D/K208G might stabilize the Ca2+ binding site which played a key role in the stabilization of PML. This study provides an efficient and user-friendly computational design scheme for the thermal stability modification of natural enzymes and lays a foundation for the modification of PML and the expansion of its industrial applications.

Molecular design and biological activity analysis of antimicrobial peptide RIKL / 生物工程学报

Natural antimicrobial peptides have strong bactericidal activities. An obstacle of the development of antimicrobial peptides resides in the difficulty of developing peptides with high biocompatibility. In this study, molecular dynamics analysis was employed to assess the structural characteristics and biological activities of peptides. A (RXKY)2(YRY)2 structure was used as a template to design an antimicrobial peptide RIKL of high-efficiency and low-toxicity, where X represents Ile and Y represents Leu. The secondary structure of the antimicrobial peptide was detected by circular dichroism (CD), and the structures of RIKL in water and in POPC/POPG membrane environment were measured using molecular dynamics. The biological activity of RIKL was further studied by assessing its antimicrobial activity, hemolytic activity, eukaryotic cytotoxicity, and salt ion stability. CD results showed that RIKL presented an α-helical structure in a simulated bacterial membrane environment. Molecular dynamics simulation predicted that the secondary structure of RIKL could be partly retained in water and POPG environment, while this secondary structure was weakened in the POPC environment. Antimicrobial test suggested that RIKL had high antimicrobial activities, and the geometric mean of the Minimum Inhibitory Concentration (MIC) was 3.1 μmol/L. The hemolysis indicated that RIKL had no hemolytic activity within the detection range, and cytotoxicity test suggested the cytotoxicity of RIKL was low. Stability test showed that RIKL maintained antimicrobial activities under different pH, serum concentrations and salt environments. Based on the above results, RIKL has high cell selectivity and has the potential as a highly effective antibacterial drug.

Identification of small molecule sirtuin-1 inhibitors for treating acute myeloid leukemia based on molecular docking, quantitative structure-activity relationships and molecular dynamics / 药学学报

The purpose of this study was to discover novel inhibitors of sirtuin-1 (SIRT1) that could be used in the treatment of acute myeloid leukemia (AML). Eight potential SIRT1 inhibitors were identified from 231 511 natural drug-like molecules by virtual screening-based molecular docking and molecular mechanics-generalized Born surface area (MM-GBSA) calculation of binding free energies. Using existing SIRT1 inhibitor molecules as training and test sets, a series of quantitative structure-activity relationship models were established, and the best quantitative structure-activity relationship (QSAR) model was used to predict the IC50 of these 8 potential inhibitor molecules for SIRT1. Subsequently, molecular dynamics simulations were performed to verify the binding mode and stability of these complexes of potential inhibitors and SIRT1 protein. Finally, the activity of these potential SIRT1 inhibitors was verified by cell proliferation assays of OCI-AML2, OCI-AML3 and MV4-11 cells and SIRT1 enzyme activity assays, and it was found that 5 compounds could inhibit AML cell proliferation. Among them, the most active compound, ZINC000001774455, had an IC50 of 2.29 ± 0.09 μmol·L-1 with OCI-AML2 cells, and at a concentration of 1 μmol·L-1, the inhibitory ratio of this compound on SIRT1 protein activity was 65.33%. ZINC000001774455 can be used as a lead compound for the development of new AML treatments.

Computation-aided design of the flexible region of zearalenone hydrolase improves its thermal stability / 生物工程学报

The zearalenone hydrolase (ZHD101) derived from Clonostachys rosea can effectively degrade the mycotoxin zearalenone (ZEN) present in grain by-products and feed. However, the low thermal stability of ZHD101 hampers its applications. High throughput screening of variants using spectrophotometer is challenging because the reaction of hydrolyzing ZEN does not change absorbance. In this study, we used ZHD101 as a model enzyme to perform computation-aided design followed by experimental verification. By comparing the molecular dynamics simulation trajectories of ZHD101 at different temperatures, 32 flexible sites were selected. 608 saturated mutations were introduced into the 32 flexible sites virtually, from which 12 virtual mutants were screened according to the position specific score and enzyme conformation free energy calculation. Three of the mutants N156F, S194T and T259F showed an increase in thermal melting temperature (ΔTm>4 °C), and their enzyme activities were similar to or even higher than that of the wild type (relative enzyme activity 95.8%, 131.6% and 169.0%, respectively). Molecular dynamics simulation analysis showed that the possible mechanisms leading to the improved thermal stability were NH-π force, salt bridge rearrangement, and hole filling on the molecular surface. The three mutants were combined iteratively, and the combination of N156F/S194T showed the highest thermal stability (ΔTm=6.7 °C). This work demonstrated the feasibility of engineering the flexible region to improve enzyme performance by combining virtual computational mutations with experimental verification.

Effect of graphene oxide on the function of erythrocytes / 生物工程学报

The biocompatibility of nanomaterials has attracted much attention. Graphene oxide (GO) is a nanomaterial widely used in biomedicine, but its toxicity can not be ignored. In this study, the effect of GO on the blood system (the hemolysis rate, the fragility of erythrocyte, and acetylcholinesterase activity) was systematically investigated. The results showed that the hemolysis rate of erythrocytes was lower than 8% when the GO concentration was below 100 μg/mL (P5 μm (LGO) increased the activity of acetylcholinesterase by 42.67% (P<0.05). Then molecular dynamics simulation was used to study how GO interacted with acetylcholinesterase and increased its activity. The results showed that GO was attached to the cell membrane, thus may provide an electronegative environment that helps the hydrolysate to detach from the active sites more quickly so as to enhance the activity of acetylcholinesterase.

Structure-based optimization and design of CRISPR protein xCas9 / 生物工程学报

Streptococcus pyogenes Cas9 (SpCas9) has become a powerful genome editing tool, but has a limited range of recognizable protospacer adjacent motifs (PAMs) and shows off-target effects. To address these issues, we present a rational approach to optimize the xCas9 mutant derived from SpCas9 by directed evolution. Firstly, energy minimization with the Rosetta program was applied to optimize the three-dimensional structure of Cas9 to obtain the lowest energy conformation. Subsequently, combinatorial mutations were designed based on the mutations sites of xCas9 acquired during the directed evolution. Finally, optimal mutants were selected from the designed mutants by free energy ranking and subjected to experimental verification. A new mutant yCas9 (262A/324R/409N/480K/543D/694L/1219T) with multiple PAM recognition ability and low off-target effects was obtained and verified by DNA cleavage experiments. This mutant recognizes the NG, GAA and GAT PAMs and shows low off-target DNA cleavage activity guided by mismatched sgRNA, thus provides a gene editing tool with potential applications in biomedical field. Furthermore, we performed molecular dynamics simulations on the structures of SpCas9, xCas9 and yCas9 to reveal the mechanisms of their PAM recognition and off-target effects. These may provide theoretical guidance for further optimization and modification of CRISPR/Cas9 proteins.